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We study the competition between interactions and disorder in two dimensions. Whereas a noninteracting 
system is always Anderson localized by disorder in two dimensions, a pure system can develop a Mott gap 
for sufficiently strong interactions. Within a simple model, with short-ranged repulsive interactions, we show 
that, even in the limit of strong interaction, the Mott gap is completely washed out by disorder for an infinite 
system for dimensions D < 2, leading to a glassy state. Moreover, the Mott insulator cannot maintain a broken 
symmetry in the presence of disorder. We then show that the probability of a nonzero gap as a function of system 
size falls onto a universal curve, reflecting the glassy dynamics. An analytic calculation is also presented in ID 
that provides further insight into the nature of slow dynamics. 

PACS numbers: 



I. INTRODUCTION 

It has been a dream in condensed matter physics to describe 
the quantum phase transition between localized and itinerant 
electrons, as it is a reflection of the basic concept of wave- 
partcle duality in a quantum many body system. Itinerancy 
mirrors the wave aspect, localization the particle aspect. In 
one-particle quantum mechanics, wave and particle descrip- 
tions are dual of each other, and there is no fundamental dis- 
tinction between them. Coherent superposition of waves are 
packets that act like lumps of energy, or particles. Yet, in a 
many particle system one believes that the metallic state, de- 
scribed by a non-normalizable wave function, is separated by 
a quantum phase transition from normalizable localized states, 
where particles are tied to spatial centers. 

Band theory proposes a sharp distinction between metals 
and insulators. Although a typical eigenstate carries current, 
the totality of electrons in a filled band cannot. Paradoxically, 
in spite of the quantum mechanically coherent and extended 
nature of each electronic eigenstate, the system is an insula- 
tor. The many particle wave function of an insulator is a Slater 
determinant of Bloch functions of a filled band. Alternately, 
the same determinant can be rewritten as a determinant of lo- 
calized Wannier functions. This is a manifestation of wave- 
particle duality. 

An interaction driven insulator, or a Mott insulator, can be 
an insulator even if the band is half-filled and can be due to 
a local repulsive, at most a few body, interactions. While 
this can lead to a collective localized state, this mechanism 
is vastly different from the non-local statistical constraint en- 
forced by the Pauli exclusion principle, as in a band insula- 
tor of non-interacting electrons. Mott insulators are similar 
to classical insulators. Without quantum mechanics, at zero 
temperature, a system of electrons will assume the configu- 
ration of the lowest potential energy due to interactions, and 
because of the harmonic restoring force they will not conduct 
in response to an applied electric field. The lowest energy 
state is likely to be a broken symmetry state with crystalline 
order. A classical insulator is a localized state stabilized by 
interactions. 

There is another remarkable alternative, the Anderson insu- 
latojP^. The non-interacting electronic eigenstates may them- 



selves localize due to a random potential, and if the Fermi 
energy is situated within the localized states, the system is 
an insulator Like a band insulator, quantum interference lo- 
calizes a particle due to interference of time reversed paths, 
another manifestation of wave-particle duality. A priori it is 
not clear when this physical situation realizes, as the role of 
interaction becomes more and more important as the system 
approaches localization. Nonetheless, we would like to show 
that in certain circumstances the opposite may be true, that is, 
disorder dominates, however weak it may be. 

One of the mechanisms by which an interacting system 
without disorder may become insulating is by opening a gap 
in the excitation spectra by breaking symmetries, such as spin 
and charge density waves that are particle-hole condensates. 
This mechanism provides a definition of Mott insulators, in 
the sense that a half-filled band could insulate. Whether or 
not all Mott insulators must be accompanied by a broken sym- 
metry has been the subject of some recent debatJ^Ei Rather 
than addressing this issue, we shall assume that there is a bro- 
ken symmetry in the Mott state, which is often the case, and 
in fact it is the reason for its existence. We emphasize that 
broken symmetry is a general concept for which correlation 
effects are sine qua non. Thus, the mechanism itself must not 
be identified with a Hartree-Fock approximation. 

The nature, and even the existence, of a 2D metal-insulator 
transition in low disorder Si-metal-oxide-semiconductor field- 
effect transistors has remained controversial since it was first 
reported^, despite considerable experimental and theoretical 
efforffl The fundamental difficulty is in understanding the 
complex interplay between strong interactions and quenched 
disorder Noninteracting electrons (or even a Fermi liquid) are 
localized by any amount of disorder in two dimensions (2D^, 
which implies that if a metallic phase is found in experiments, 
it must reflect a non-Fermi liquid ^. 

Due to their complexity, a principled analysis of systems 
involving both strong interactions and disorder is necessary 
to understand what sorts of qualitative behaviors may result 
from these two basic ingredients. In this work, we will pro- 
vide such an analysis, albeit in a simple model. The paper is 
organized as follows. In Section II, we introduce the model 
and show, quite rigorously, that this model does not have a 
true metal-insulator transition in the presence of disorder The 
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analysis is performed in the strong interaction limit where the 
Mott gap is the largest and should be the most resistant to 
the onslaught of disorder. In Section III, we study, through 
numerically constructed ground states, the probability for a fi- 
nite system to have a non-zero gap, and show that this quantity 
falls onto a universal curve. From this result, we then show 
that the system possesses glassy quantum dynamics. In Sec- 
tion IV, we then specialize to the case of one dimension and 
study analytically the disorder averaged ground state density 
as a function of chemical potential. We compare this to the 
density achieved from a rapid temperature quench to illustrate 
anomalously slow dynamics even in II?. Finally, we conclude 
in Section V with a brief discussion of the implications of our 
results. 



II. THE MODEL 



variables with zero mean and variance cP' . Consider the case 
where we drop the hopping terms and we are left simply with 
a classical model since the A term clearly commutes with H^. 
We do this not as an approximation, but to prove the point that 
even in the limit that the Mott insulator has the best chance of 
surviving, any amount of randomness destroys the Mott ga£. 
We are strongly motivated by an argument offered for IE™. 
As A is lowered, the kinetic energy will cause the walls to 
fluctuate, but this won't change the fact that the gap has been 
destroyed by randomness, nor will it restore the symmetry. As 
long as the disorder remains finite, it is difficult to believe that 
the system will ever reach a true metallic state. 

Using rii — Sz{i) + 1/2, we see that we have, in fact, an 
Ising model. Redefining Sz{i) —Sz{i) at every other site, 
and remembering that the Vi are symmetrically distributed 
with zero mean, gives simply the classical ferromagnetic ran- 
dom field-Ising model (RFIM): 



Consider the simplest 2D fermion model that has a bro- 
ken symmetry in the insulating state: spinless fermions on a 
square lattice with the Hamiltonian 
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The sums are over nearest neighbors of a bipartite lattice 
with sublattices A and B. The symmetry under the opera- 
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in A, while Ci c| in B, ensures half-filling, 
with rii — clci being the density operator at site i. This model 
is well studied^ in ID where the system is a Mott localized 
insulator in our sense, a charge density wave state, with a gap, 
g, for A > 1; 5 — > A, for A 3> 1. For A < 1, the system is 
metallic, therefore A 1 is the location of a metal-insulator 
transition. The same transition must obviously be present 
in 2D on general grounds, as the fluctuations are weaker in 
higher dimensions. 



FIG. 1: Example of a domain wall separating a region of down spins 
(above) from up spins (below). The dotted line is an example of a 
(gapless) comer excitation in which no extra bonds are broken. 

Addition of a random potential gives H = Hp + H,. with 
Hr = — Vini - where Vi are independent Gaussian random 
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Although 2D is the marginal dimension for the Imry-M^^l 
argument, it has been rigorously shown by Aizenmann and 
that there is no long-range order for arbitrarily weak 
randomness. Moreover, since the ground state will generically 
be a disordered spin configuration dependent upon the partic- 
ular realization of the random potential, the excitations of the 
system result from moving domain walls between the up and 
down spin regions. As a result of the disorder, the domain 
walls will in general be rough, and the elementary excitations 
then consist of moving corners because no extra bonds are 
broken (see Fig. [T]i. The energy cost of moving comers de- 
pends only on the random field configuration and not on A. 
Thus, for an infinite system, the excitation spectrum will be 
essentially gapless for a continuous distribution of Vi, which 
we have assumed. Note that the corner excitations are not nec- 
essarily the first excited states, but they are sufficient to prove 
the existence of gapless excitations. 



III. TWO-DIMENSIONAL GROUND STATES AND 
GLASSY DYNAMICS 

Having argued against the existence of a metal-insulator 
transition in the presence of disorder, we turn to an analy- 
sis of the ground state in 2D and provide evidence for glassy 
quantum dynamics. As already mentioned, D = 2 is special 
for the RFIM. In \D, the Imry-Ma argument easily gives that 

ordered domains have size Lc ^ (^) . But since D = 2 
is the marginal dimension, no information can be gleaned 
from the simple Imry-Ma argument. Rather, it is necessary 
to study the energy gained upon allowing the domain walls 
to roughen. Such a calculation was performed by Binder'^ 

with the result,tSIBl that Lc ^ exp yl (^)^ , where A is a 

constant. Thus, domains in two dimensions are exponentially 
larger than those in ID. It is important to note that this equa- 
tion relates the typical size of domains to the ratio of exchange 
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to random field energies. 

This leads us to ask the following question: Given a finite 
system of size L x L, what is the probability it possesses an 
energy gap? For small L ^ L^., a non-zero gap would be 
fairly likely, but as L increases towards the crossover length, 
Lc, given above, the probability of finding a gap should de- 
crease because the system may now contain multiple ordered 
domains separated by rough walls. 



A. Numerically Computed Ground States 

We answer this question numerically by computing exact 
ground state configurations of the RFIM for different system 
sizes and disorder strengths. The general method'^ for find- 
ing ground states of random field systems (or even random 
bond systems without frustration), is based on a mapping to 
an equivalent minimum cut network flow problem. In a net- 
work, nodes (i.e. lattice sites) are connected by directed links 
with finite capacity, signifying the maximum possible flow be- 
tween neighboring nodes. Two additional sites (dubbed the 
source and sink) are augmented to the lattice, and each site 
of the RFIM is connected by a directed link to one of the 
two external sites, depending on the sign and strength of its 
random field. The flow capacity between neighboring lattice 
sites is determined by the exchange energy. Dividing the net- 
work in two, with source and sink on opposite sides of the 
division, defines a cut. With this construction, the minimum 
of the capacity across all possible cuts provides the ground 
state energy of the RFIM, while the minimum cut itself deter- 
mines the spin configuration. Using the equivalence between 
the minimum cut capacity and the maximum flow through the 
network, known as the max flow-min cut theorem'^, allows 
one to simply calculate maximum flows. To do this, we em- 
ployed the efficient push-relabel codi^ which enables us to 
get good statistics for moderate system sizes. To illustrate 
typical ground state domain structures for various disorder 
strengths, we computed minimum cuts explicitly through the 
Edmonds-Karp algorithm'^, which are displayed in Fig.pl 

For each chosen value of cr/A and InL, we calculate the 
fraction of realizations of the Vi that contain a domain wall 
in their ground state. To detect a domain wall, the exact 
ground state energy computed via the network flow model is 
compared with the minimum energy of the two ferromagnetic 
states (all up or all down). If the ground state energy is lower, 
there must be a domain wall. If not, and the ground state en- 
ergy equals the lower energy ferromagnetic state, then there 
must not be a domain wall. This procedure avoids having to 
examine the spin configuration explicitly. In other words, we 
will assume that the presence of a domain wall implies gap- 
less excitations even for a finite system. On the other hand, if 
the ground state is purely ferromagnetic, the excitation energy 
will be non-zero and of order A. Figure [3] shows our results. 
As expected, the probability of domain wall formation grows 
upon increasing both In L and cr/A. The surprising feature is 

that when we define the x-axis to be a; = InL — A (^)^, 
all points collapse onto a single universal curve. We find 
A = 1.8, in good agreement with a previous study which re- 




FIG. 2: Typical ground state configurations for different disorder 
strengths and system sizes. Top: L = 30 x 30 and A/ a = 3. 
Bottom: i = 20 x 20 and A/a = .95. 

ported ^ = 2.1 ± .plat flie special value P{g = 0) = 1/2. 
The collapse onto a single curve enables us to gain informa- 
tion about the slow transition to a disordered ground state even 
for values of a/ A for which the relevant system sizes are far 
too large to be studied numerically. Most notably, not only 
does the typical size for observing a disordered ground state 
(i.e. the L for which P{q = 0) = 0.5) scale as predicted 
iiP, the entire distribution scales in precisely the same way. 
This is surprising because it might have been expected that 
systems with particularly weak randomness cross over to dis- 
ordered ground states more 'slowly', i.e. over a much broader 
range of InL, in addition to reaching P{g = 0) = 0.5 at a 
larger length scale. Note also that the gap probability is also 
not symmetric about P{g = 0) = 0.5. 

The crossover from a generically ordered {P{g = 0) ^ 0) 
to disordered {P{g = 0) ~ 1) ground state occurs quite 
slowly, over nearly two decades. The large range of L for 
which samples are neither generically ordered nor disordered 
indicates the possibility of large sample to sample fluctua- 
tions. Some systems may have remnants of a Mott gap, while 
others are localized by disorder In addition, the large L 
regime will likely possess many low-lying energy states, each 
requiring the reorganization of large numbers electrons, im- 
plying significant metastability and glassy dynamics, akin to 
experiments on 2D low-mobility Si inversion layers^' ^^^^ 
although the interactions considered here are short-ranged. 
Such glassy behavior has also been found theoretically in a 
similar model of spinless fermions on a Bethe lattice^SI It is 
quite remarkable that contrary to expectations the system be- 
haves more like an "Anderson insulator" despite strong inter- 
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FIG. 3: (Color online) Probability to find an energetically favorable 
domain wall in a finite system of linear size L. The values of A/ct 
range from 0.2 to 1, while the values of L are 20, 40, 80, 160, and 
320. Symbols represent different fixed system sizes with varying 
disorder strengths. For each choice of A /o and L, we averaged over 
3,000 realizations of disorder. We find that A = 1.8 provides the 
best data collapse. The black line is the fit to an asymmetric sigmoid 
(see text). 

actions. 



B. Glassy Dynamics 

The signature of glassy quantum dynamics lies in the long 
time-scale for relaxation to the ground state. We can, in fact, 
learn something about the dynamics by interpreting Figure 3. 
The distribution P{g = 0) can be regarded as the cumula- 
tive distribution function (cdf) to find an ordered domain of 
characteristic size L or smaller, since this distribution func- 
tion reflects the existence of a domain wall up to the scale L. 
The curve fits an asymmetric sigmoid f{x) of the form 

/(^) = 7 / ,,,,e ^ (3) 

where x = \rL{Le~^^^/'^^ ). The best fit to the data shown 
in Fig. |3]is 9 = 0.31, xq = 1.37, and ^ = 0.29. Tak- 
ing a derivative to find the L distribution, P{L), results in 
P{L) ~ L^^l^ for large L. If we define the imaginary part of 
the frequency dependent local susceptibility corresponding to 
the density riiit) (Heisenberg operator) to be x"('-^)' then 

X"H ~ y P(L) 8{uj - cc-oe-^^"). (4) 

The (5-function signifies that at a frequency the quantum 
tunneling rate corresponding to that frequency is sampled by 
a cluster of size L, where the exponent a requires micro- 
scopic calculation and is left undetermined in the present phe- 
nomenological analysis. The quantity ujq is the attempt fre- 
quency in the many dimensional WKB theory. Thus, it is easy 



to show that as 0, 

The exponent -0 = 1 + (V^ ^ > 1' long as a 7^ 0; 
a = is highly unlikely because that would imply that the 
action corresponding to the tunneling rate is independent of 
the size of the cluster, L. It should be interesting to check ex- 
perimentally that the noise power spectrum does follow this 
1 /(jj-Iaw with a logarithmic correction, signifying glassy dy- 
namics. 



IV. ANALYTIC RESULTS IN ONE DIMENSION 

It is instructive to consider the same model of spinless 
fermions in one dimension"' where it is possible to compute 
disorder averages of thermodynamic quantities analytically. 
Recall that in \D, the Imry-Ma argument gives that the or- 
dered domains have characteristic size L^. ^ (^) , so the 
ground state is "disordered" and heterogeneous. To illustrate 
the origin of glassy dynamics present even in 1_D, we will 
calculate both the disorder averaged ground state density pQ 
as well a quenched density pQ that would be obtained from 
an infinitely fast temperature quench. This quenched state is 
obtained through the sequential filling of the lowest available 
energy levels up to the chemical potential p.. We will be inter- 
ested in how the density profiles vary with p, so we must relax 
the condition of half-filling. In addition, we will work in the 
limit of large nearest-neighbor repulsion, A 00, so that the 
particles must be separated by at least one empty lattice site. 
In other words, the fermions may be regarded as hard-core 
dimers (see Fig. |4]l. We emphasize that these dimers should 
not be confused with valence bonds connecting two neighbor- 
ing sites, for which the word "dimer" is frequently also used. 
The hard-core dimer constraint, along with the random on-site 
energies, induces a geometrical frustration between compet- 
ing particle configurations. As a result, the ground state den- 
sity profile is a complex structure that incorporates the pref- 
erence for low energies while still respecting the hard-core 
constraint. 



A. Heuristic Analysis of Dimer Frustration 

Before calculating pQ and pq, we give a simple argument 
to show that these two quantities differ by a finite amount 
for all pg > 0. In what follows, the on-site energies will 
be drawn from a uniform (i.e. rectangular) distribution, de- 
noted by i?(e), between zero and one. The uniform distribu- 
tion will simplify the analytics and preserve the consistency 
of the dimer (i.e. large A) limit. When p=Q, all sites have 
positive energy, so the ground state is an empty lattice. When 
p = 1, all sites are attractive, but the ground state has a com- 
plicated structure due to the competing effects of the random 
on-site energies and the hard-core constraint. Therefore, we 
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FIG. 4: (Color online) Schematic of a scenario where the quenched 
configuration occupies site k, hence blocking sites fc— 1 and fc+1, but 
the ground state foregoes occupancy of site k in favor of the flanking 
sites if El + E2 < E. The upper (blue) configuration is the ground 
state and the lower (red) is the quenched. 



B. Exact Solution for Ground State Density 

We now derive the disorder-averaged ground state density 
Pa for a ID lattice of dimers at chemical potential and with 
on-site energies uniformly distributed between zero and one. 
The analysis follows that of Fonk and HilhorsP who con- 
sidered the problem with a different energy distribution not 
easily generalizable to a non-zero chemical potential. Define 
EI{E1) to be the ground state energy of the first k sites sub- 
ject to the constraint that a particle is present (absent) at site 
k. These quantities can easily be seen to obey the recursion 
relations, 

El = e, + El_, (7) 
El^Tnin{El_,,El_,) (8) 

where Ck is the random energy at site k. Defining difference 
variables = El — E^ and subtracting the above equations 
gives a simple recursion relation for ^^j, 



will focus on /i between zero and one. In the following ar- 
gument, we will absorb /i into the on-site energies which will 
instead be uniformly distributed between — /i and 1 — fj,. 

Consider a finite lattice of L sites. Let the energy at site k 
be i? < and the energies at sites fc — 1 and fc + 1 be i?i and 
E2, respectively. We calculate the probability, ^'syjyjtch' 
in the ground state, placement of a particle at k is forfeited in 
favor of the occupation of sites fc — 1 and fc + 1, despite site fc 
having the lowest energy. The scenario is depicted in Figure 
|4] This is the simplest way geometric frustration may cause 
the ground and quenched states to differ because poik) — 
but pqik) = 1. Averaging over the value of E, 



switch 



P{Ei +E2< E)dE 



(6) 



Since the sites fc — 1 and fc+1 must also be attractive, 
P{Ei + E2 < E) ^ J^dEiJ^dE2e{E~{E,+E2)) 
which equals £'^/2. Plugging this into (j6|l then gives that 
-Pswitch = Multiplying Pswitch by the average num- 

ber of sites with E < 0, i.e. pL, gives the average number of 
these switches in a finite system of size L. Setting this equal 
to unity gives the chemical potential at which we expect the 
first switch: pc ^ L~^/^. Clearly, as L —> 00, /ic ^ 0. Also, 
since each such switch increases the ground state density rel- 
ative to the quenched, we also have that 5 = pq — pq ^ p^ 
for small /i, so (5 = only at p = 0. The above argument ne- 
glects the contributions of sites fc — 2 and fc + 2 etc. but these 
effects are higher order in p and hence can be neglected for 
/Lt <C 1. We will see this behavior of 5 reproduced precisely 
from the exact result. 



Cfe = efc + min {0,-£_k-i) 



(9) 



By averaging over the on-site energy distribution R{e), the 
recursion relation is readily transformed into an integral recur- 
sion relation for P(^), the distribution function of ^, 
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(10) 



The fixed point distribution of P{£,) will contain the re- 
quired information for bulk quantities, so we can drop the sub- 
scripts on P(^). With i?(e) = 0{e + p)9{l - p - e), we see 
that P{£,) — for ^ > 1 — /i and hence also for ^ < — 1. Then 
there are three distinct regions to consider: 

Region 1: — 1 < ^ < —p 



P(0 = PiOde 



(11) 



Region 2: < ^ < 



P(0 = y° P(C')d$' + y ' ' PiO^e - 1 (12) 



Region 3: < ( < I - p 



P{e)de (13) 



If we know the solution in region 3, we can integrate to find 
the solution in region 1. Region 2 has a flat value of 1 (since 
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P{£,) is a normalized probability distribution). We convert 
( fTSj l to the differential equation 



dm 



(14) 



which can be reduced to two coupled ODEs by the replace- 
ment Q{£,) = P{1 — /i — ^). The resulting solution for P{^) 
in region 3 is 



F(0=cose 



sin (V)- cos (V) 
sin(V)+cos(V) 



sinC (15) 



There is an undetermined multiplicative constant fixed by 
requiring that P{^) integrates to 1 . From the form of the equa- 
tion in region 3, we see that this is equivalent to requiring that 
P(l — /i) + Jq ^ P(^')d^' = 1. The necessary constant turns 
out to be unity. One can then directly integrate to find the 
solution in region 1 : 



P{0 



2 sin 



sm 



cos 



S±i£ 
2 



> 2 ) 



(16) 



Eqs. ( 15 I and ( 16 1, along with P{Cj = 1 in region 2, com- 
prise the required solution of the integral equation. 

We now use the derived form of P(^) to solve for the 
disorder-averaged ground state density. This can be found, 
again following^*, by defining E^{E^) to be the minimum 
energy of the entire system (not just the left half), subject to 
the constraint that a particle is absent (present) at some site 
k deep in the bulk. The average density will be given by 
\ — P{E^ < E^). A similar recursive calculation, this time 
including sites to the left and right, leads to 



p(£;o < pi) 



' ''deiP(ei) ""d^Pfe) X (17) 

-1 J -I 



j ''de0(-e-min[O,-ei] -min[0,-$2]) 
The theta function can be split up into four cases 



de9{ 



e — min[0, — ^i] — min[0, —£,2]) 



(18) 



0(^1)0(6) min(l, ^1 + ^2 + A*) + 

and each term integrated with the form of P(^) derived above. 
After some lengthy but straightforward algebra, and using 
Pq = 1 — P{E^ < pi), we find the remarkably simple result 



1 



PG = 



(19) 



1 + CSC fj. 

It is important to remember that this result is valid in the 
regime < /i < 1. Other quantities such as the average 
energy per site can be calculated from P(^), should they be of 
interest. 
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FIG. 5: (Color online) (A) Plot of the analytic forms of the dis- 
order averaged ground (red solid line) and quenched (blue dashed 
line) state densities. Both initially rise linearly from zero, but the 
quenched density peels off for higher /i. Note that pa is bounded 
below by pQ and the two functions only intersect at zero. (B) Plot of 
S ^ PG~ PQ vs. PC- 



C. Quenched State Density 



We now turn to the calculation of the quenched state den- 
sity. To do this, we will use the formalism of random sequen- 
tial adsorption (RSA).'^In particular, we will use the dynamic 
formulation of RSA with a random distribution of binary ad- 
sorption rates^** in the determination of the quenched density. 
The reason is that, since the random energies on distinct sites 
are uncorrected, the process of sequentially filling the deep- 
est energy minima is identical to an RSA process^^. However, 
with /i < 1, some sites are repulsive and hence have an "on 
rate" of zero. Therefore, we need to consider an RSA pro- 
cess with two adsorption rates, a and /3, take the limit f3 ^ 
(while a remains arbitrary), and look for the < 00 density. 

The adsorption rate of site n will be denoted by a„ and 
the probability that site n is occupied by a particle at time t 
is Pn{t)- This probability should be thought of as an average 
over different realizations of the adsorption process for a fixed 
choice of the a„. Using estabUshed formalisrrP^, p„ varies in 
time as 
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dpnjt) 

dt 



an exp(-a„t)Q„_^i(3+_i 



(20) 



where the Q's are time-dependent and obey 



dQn 

dt 



-an exp{-ant)Qn+i 



(21) 



with a similar equation for (5+ except with the replacement 
n+1 — > 1. Since in our calculation of the ground state den- 
sity, we have chosen the on-site energies to be uniformly dis- 
tributed between zero and one, the chemical potential ^ gives 
the fraction of attractive sites that therefore have adsorption 
rate a 0. Thus we choose the q;„'s to be random variables 
equal to a with probability /i and equal to zero with probabil- 
ity 1 — /i. The rest of the sites are repulsive and hence have an 
adsorption rate of zero. The sequential filling of lowest energy 
minima will then be mimicked by the dynamic RSA process. 

To find the average density of the quenched configuration, 
Pq, we need to average over the q;„'s and take the t — > oo 
limit. Since Q'+i only depends on sites m > n + 1 (and 
Qn-i only on TO < n — 1), the average over q;„, denoted by 
(.), simply factorizeJ^ 



dt 
d{Q) ^ 
dt 



— aexp(— ai) (Q) 
— /iacxp(— ai) {Q) 



(22) 
(23) 



where we have used (a„ exp(— Q!„f)) = /iQ;exp(— at) and 
the fact that the Q's become independent of position after 
averaging, due to translational invariance. Clearly, we then 



have {Q) = cxp [p{e — 1)] and upon integration of (22 1 



fori 



we find that the quenched density is given by 



(24) 



Pq rises linearly from zero at /i and saturates at /i = 1 to 
the "jamming" density of dimers ~ .432. 



D. Discussion 



predicted earlier in our heuristic argument. We also plot 8 vs. 
Pg in Figure |5] to show that 6 remains quite small until the 
lattice reaches quarter filling {pQ — .25). 

The picture that emerges from our exact solution leads to 
a few general conclusions. First, the quenched and ground 
states of hard core particles in disordered landscapes always 
differ in an extensive fashion, except at = 0. As /i is 
increased, the ground state evolves through local rearrange- 
ments describable as "micro" first-order transitions^**. Despite 
this, there are two qualitatively different regimes where 5 can 
be either large or small. In the small 5 regime, the ground state 
is kinetically accessible, whereas the large 5 regime is char- 
acterized by ubiquitous metastability. If the system becomes 
stuck in one of these metastable configurations, the geomet- 
ric frustration will cause anomalously slow evolution towards 
the ground state. From Fig|5] it is clear that the symptoms of 
glassy dynamics become more pronounced once the system 
reaches quarter filUng. 



V. CONCLUSION 

Having established the main points of the paper it is use- 
ful to make a few educated guesses that should be of interest 
to future work. We reiterate that as long as the disorder re- 
mains finite, it is unlikely that the system will reach a true 
metallic state. If the picture described here is generic (ignor- 
ing long range Coulomb interaction), there would not be a 
true metal-insulator transition in 2D. However, there should 
be a crossover scale below which the system will appear to 
have an insulating Mott gap. In contrast, the corresponding 
?>D case could exhibit a genuine quantum phase transition be- 
cause the Imry-Ma argument leaves open the possibility of a 
broken symmetry state signifying a Mott insulator, and the 
argument for the proliferation of low energy excitations in 
D < 2 cannot apply. There are many systems where patchy 
gapped states or a filled Mott gap appear to be important as 
in underdoped high temperature superconductors or frustrated 
magnets. It should be interesting to examine the role of dis- 
order from the present perspective. Since we have provided 
a relatively complete characterization of the ID model in the 
strong interaction limit, it would also be interesting to per- 
form a DMRG simulation where one can explicitly tune the 
interaction strength. 



The two densities ( 19 1 and (24 1 are plotted in Figure |5] and 
their difference 6 is plotted below, pg provides a lower limit 
for Pg and both rise linearly from zero for small values of 
p. For larger p, pq peels off due to the jamming caused by 
irreversible adsorption. When we expand & — Pc ^ Pq for 
small /i, the first non-zero term is /i'^/6, identical to what was 
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